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We study the spin-crossover molecule Fe(phen) 2 (NCS )2 using density functional theory (DFT) 
plus dynamical mean-field theory, which allows access to observables not attainable with tradi¬ 
tional quantum chemical or electronic structure methods. The temperature dependent magnetic 
susceptibility, electron addition and removal spectra, and total energies are calculated and com¬ 
pared to experiment. We demonstrate that the proper quantitative energy difference between the 
high-spin and low-spin state, as well as reasonably accurate values of the magnetic susceptibility 
can be obtained when using realistic interaction parameters. Comparisons to DFT and DFT-I-U 
calculations demonstrate that dynamical correlations are critical to the energetics of the low-spin 
state. Additionally, we elucidate the differences between DFT+U and spin density functional the¬ 
ory (SDFT) plus U methodologies, demonstrating that DFT-I-U can recover SDFT-I-U results for 
an appropriately chosen on-site exchange interaction. 


The combination of density functional theory (DFT) 
and dynamical mean-field theory (DMFT) is now estab¬ 
lished in condensed matter physics as a successful theory 
of materials with strong local electron correlations.^’^ Ini¬ 
tially devised as a theory of extended (infinite) systems, 
the method has been extended to finite systems^^® and 
has been used to demonstrate that many-body effects 
are important for ligand binding on the active center of 
protein myoglobin and haemoglobin.As compared 
to traditional highly accurate quantum chemical meth¬ 
ods DFT-I-DMFT provides many advantages including 
excited-state properties, nonzero temperatures and treat¬ 
ment of arbitrary strength of local correlations. Also, 
because the computational cost scales linearly with the 
number of symmetry-inequivalent correlated atoms, the 
method can be used to treat molecules containing many 
transition metal or actinide atoms. However, its broad 
applicability and quantitative effectiveness in the quan¬ 
tum chemical context is not yet fully established. 

Here we apply the DFT-I-DMFT method to study spin- 
crossover complexes: molecular species that change spin 
state upon increase of temperature or other changes in 
environment. Spin crossover molecules provide an impor¬ 
tant challenge to theory, requiring both accurate energet¬ 
ics and the ability to treat excitations in a situation that 
(because of the spin) necessarily involves strong electron 
correlations. Insights gained from study of spin crossover 
materials could potentially be useful for the design of thin 
films^^’^^ or single-molecule^^ spintronic devices. 

This paper provides a comprehensive DFT-I-DMFT de¬ 
scription of Fe(phen) 2 (NCS) 2 ,^® a member of an exten¬ 
sively studied and still expanding family of spin-crossover 
complexes based on Fe(H).^®’^^ We compute the mag¬ 
netic susceptibility, electron addition and removal spec¬ 
tra, and total energy, finding results in good agreement 
with experimental data when using realistic interaction 
parameters. Our analysis enables us to infer that the 
metal-to-ligand bond length is the control parameter 


of spin transition. We explore the sensitivity of vari¬ 
ous observables to the double-counting correction and 
the on-site interactions U and J, demonstrating that 
Fe(phen) 2 (NCS )2 is a useful testbed for current and fu¬ 
ture first-principles methods. Comparison of our results 
to those obtained with the Hartree approximation (ie. 
DFT-|-U^®’^®), demonstrates the importance of dynam¬ 
ical fluctuations in capturing the physics of strong hy¬ 
bridization that is present in the LS state. 

The key feature of Fe(phen) 2 (NCS )2 is the octahe- 
drally coordinated Fe(H) ion.^®’^^ As the temperature 
is increased above T* = 176A", there is an abrupt in¬ 
crease in magnetic susceptibility which is believed to 
be related to a change in the electronic configuration 
from the nominal low-spin (LS) state with no un¬ 

paired electrons to the nominal high spin (HS) state 
with 4 unpaired electrons. The average Fe-to-N distance, 
d{Fe — N), is also longer by about 0.2A in the HS than 
in the LS state.Experimental measurements estimate 
that the energy splitting between LS and HS states is 
Ehs — Eps ~ O.lSeV.^^ This difference is much greater 
than the contribution ksT ■ ln5 ~ 0.025eF to the free 
energy from the change in electronic entropy. 

Obtaining a LS-HS energy difference of the correct or¬ 
der of magnitude has proven challenging for theory. Spin 
density functional theory in the generalized gradient ap¬ 
proximation (GGA) level overestimates the stability of 
low-spin state while Hartree-Fock theory incorrectly pre¬ 
dicts HS as the ground state.^^ These considerations mo¬ 
tivated people to consider hybrid functionals, which in¬ 
terpolate between DFT and HF energies and therefore 
can be tuned to obtain the desired energy difference. 
However, the the amount of exact exchange was found to 
be less than is normally considered reasonable.An 
extensive study of spin-crossover molecules using mod¬ 
ern density functional theory, including meta-GGA and 
hybrid meta-GGA and double-hybrid functional, shows, 
generally speaking, relative energies of spin multiplici- 
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ties are still challenging for DFT methods.The spin 
density functional plus U (SDFT+U)^®’^^ method was 
also applied to this system, but again obtaining the cor¬ 
rect energy splitting required choosing U ~ 2.5ey, much 
smaller than is believed to be relevant for Fe.^^’^® Dif¬ 
fusion Quantum Monte-Carlo method has also been ap¬ 
plied to charged spin-crossover molecules, although di¬ 
rect comparisons to experiments are not currently fea¬ 
sible for charged systems.Very recently a detailed 
quantum chemical calculation based on CASPT2 meth¬ 
ods reported an energy splitting of 0.17eV,^® indicating 
the importance of correlations in this system. 

Here we perform fully charge self-consistent 
DFT-I-DMFT calculations of the total energy using 
the method described in Ref. 30 and 31. The DFT part 
of our calculations use the Vienna ab initio simulation 
package (VASP),^^’^^ with the Perdue-Burke-Ernzerhof 
exchange-correlation functional,an energy cutoff of 
400eV and a supercell of edge length ISA. Maximally- 
localized Wannier functions (MLWF),^® constructed 
using an energy window of 26.5eV with 4.4eV of empty 
states (120 states total), were used to represent the 
so-called hybridization window and to construct the 
correlated subspace^®’^^ in which DMFT is performed. 
The hybridization window includes all 85 occupied states 
as well as the anti-bonding tt orbitals on phenanthroline 
and thiocyanate groups that hybridize with the t 2 g 
orbitals of Fe. The correlated subspace is chosen as the 
five d-like iron-centered orbitals. Following common 
practice, for the intra-d interaction we take the density- 
density part of the full on-site Coulomb interaction, 
parametrized by two independent interaction constants 
denoted as U and J (we use the form stated in Ref. 36) 
with (unless otherwise specified) U=5.0 eV and J=0.85 
eV, consistent with estimates in the literature for Fe in 
various compounds.A double counting term is 
needed to remove the on-site d interactions present in 
the Hartree and exchance-correlation functionals. We 
use the spin-independent form of Anisimov^® 

Vdc = U^{Nd-l)-J{^-l) ( 1 ) 

Park et al,^°’^^ noted that one should allow for the pos¬ 
sibility that the coefficient [/' in Eq. 1 differs from the 
coefficient U in the interaction. However in this paper we 
set U = U' everywhere except in the discussion of Eig. 4. 

The impurity model is solved using the Continuous- 
Time Quantum Monte Carlo (CTQMC) in its hybridiza¬ 
tion expansion (CT-HYB) form^^do implemented by Gull 
et al"^^ in the ALPS package.We also solved the impu¬ 
rity problem using a Hartree approximation to the inter¬ 
action. This is the DET-I-U approximation, but imple¬ 
mented using the same correlated subspace and double 
counting as in the DFT-fDMFT calculation, enabling an 
unambiguous comparison of the results obtained from the 
two methods. In both cases the whole DFT-bDMFT loop 
is iterated until the total energy difference between con¬ 
secutive updates of charge density is less than 5 meV. 



FIG. 1. Dependence of total energy on average Fe-to-N bond 
length. The DFT-I-U and DFT-I-DMFT calculations used 
U = U’ = 5.0eV and J = 0.85eV and the DFT+DMFT cal¬ 
culations were done at 387K. The two dashed lines indicate 
the experimentally measured average Fe-to-N bond length for 
LS and HS states. 


Spectral functions are obtained via analytic continuation 
of the computed imaginary time Green’s function using 
the maximum entropy method^^ and the Fe contribution 
to the magnetic susceptibility is calculated from the im¬ 
purity model spin-spin correlation functions, which are 
measured in our CTQMC calculations, as 

Xioc = {gfJ-B f [ dr(5'^(0)5'^(T)) (2) 

Jo 

where Sz is the z component of spin on the impurity site. 
The total energy (see Refs. 30 and 31 for details) is 

E*°*[p,Gioc] = 


Here is the density functional theory approxima¬ 

tion to the total energy, E^^ [p, Gjoc] is a correction to 
the DFT energy arising from the difference between the 
DFT and DFT-I-DMFT density matrices, and the inter¬ 
action energy term EP°* [Gioc] is calculated from the fre¬ 
quency dependent self-energy E and Green’s function as 


^pot = \E{iu;r.)Gioc{iu;n) 


( 4 ) 


To perform structural relaxations within our many- 
body DFT-f-DMFT theory, we first define reference struc¬ 
tures using structurally relaxed DFT calculations. The 
metastable HS (LS) state was obtained by using an ini¬ 
tialization of the Fe magnetic moment at the nominal 
high-spin (low-spin) value. We then construct a path 
between the two structures by linearly interpolating all 
atomic positions between the values found for the LS 
and HS structures and minimize the DFT-I-DMFT en¬ 
ergy along this path. We parametrize the path in struc¬ 
ture space by the Fe-N bond length. 
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TABLE I. Energy splitting and optimized average Fe-to-N 
bond length from DFT, DFT+U and DFT+DMFT calcula¬ 
tions. 


Method 

DFT DFT+U DFT+DMFT 

Expt2i’44 

Ehs — Els (eV) 

0.70 

-0.38 

0.10 

0.13 

LS d{Fe - A) A 

1.94 

1.94 

1.94 

1.97 

HS d{Ee - A) A 

2.15 

2.20 

2.18 

2.17 




FIG. 2. (a) Fe L/// edge X-ray absorption spectroscopy of LS 
and HS measured at 17K and 298K respectively from Ref. 45. 
(b) Spectral functions of empty Fe d state from analytic con¬ 
tinuation via maximum entropy method. Absolute positions 
of spectral functions are shifted to match experiment. 


Fig. 1 shows the structure dependence of the 
DFT-I-DMFT energies along with those obtained by den¬ 
sity functional and DFT-I-U methods. All three methods 
yield two locally stable structures, one with a shorter 
d{Fe — N), which will be seen to correspond to the LS 
state, and one with a longer d{Fe — N), which will be 
seen to be the HS state. 

The bond lengths and LS-HS energy splittings com¬ 
puted for the locally stable structures are given in Ta¬ 
ble I. While all methods give bond lengths in reason¬ 
able agreement with experiment, both DFT and DFT-I-U 
methods give an inadequate account of the energy dif¬ 
ferences between the LS and HS structures. DFT pre¬ 
dicts that the LS state is much too stable while DFT-bU, 
with a physically reasonable U and J, incorrectly predicts 
the HS state to be the ground state. Similar to hybrid 
functional calculations^^ and previous spin density func¬ 
tional -bU calculations,^®’^® it is possible to tune U to 
a value that reproduces the observed energy difference 
within DFT-bU, but the required U ~ 2.5 eV is unphys- 
ically small. By contrast, the DFT-bDMFT calculations 
produce a result in good agreement with experiment with 
physically reasonable interaction parameters. 

The respective electronic states of the short-bond and 
long-bond structures are found to be locally stable up to 
the highest temperatures studied ~ 1200K-, suggesting, 
in agreement with deductions from experimental data,^® 
that the electronic entropy and energetics are not enough 
to drive the observed transition. We infer from these re¬ 



FIG. 3. Temperature dependence of product yT of suscep¬ 
tibility (computed from Eq. 2 with g = 2) and temperature 
from DFT+DMFT calculations for high and low spin struc¬ 
tures at U = 5.0eV and J = 0.85eU compared to experimen¬ 
tal data.^^ 


suits that the Fe-to-N bond length is the critical vari¬ 
able; indeed calculations (see Fig. S2 in Supplemental 
Material) show that the LS to HS transition occurs when 
d[Fe — N) crosses a critical value approximately 2.10A. 
Phonon free energy will determine the actual transition 
temperature. 

From our DFT+DMFT calculations we obtain the 
many-body density matrix describing the probability of 
different configurations of the d-orbitals (see Supplemen¬ 
tal Material). We hnd, as expected, that the dominant 
configuration in the LS state has zero total spin and 
is described by an almost complete occupancy of the 
t 2 g symmetry d-states. The key issue for the energet¬ 
ics of the LS state is the correct treatment of the vir¬ 
tual charge fluctuations into the eg states, in light of 
the strong Coulomb repulsion associated with multiple 
occupancy of the Cg states. DFT predicts an Cg occu¬ 
pancy of 1.25 electrons and a relatively large hybridiza¬ 
tion energy gain. Both the tg occupancy and the hy¬ 
bridization energy gain are likely excessive due to the 
inadequate treatment of correlations in current DFT im¬ 
plementations. Alternatively, DFT+U, which adds an 
extra Hartree term, overestimates the correlation energy, 
providing an eg occpuancy of 0.99 electrons (likely too 
small) and an underestimate of the hybridization en¬ 
ergy. DMFT treats the hybridization more correctly, giv¬ 
ing an Cg occupancy of 1.15 and a reasonable value for 
the energy of the LS state. The improved properties of 
DFT+DMFT result from a proper characterization of the 
multiconfigurational character of LS state, as also found 
in quantum chemistry calculations.^^ Turning now to the 
high spin (HS) state, we find that the HS state is found 
to be the d® maximal spin configuration, with only small 
quantum fluctuations towards d® and lower total spin, 
leading to a mean d-occupancy ~ 5.3. By neglecting the 
correlated nature of the virtual hopping into the d® con¬ 
figuration DFT+U allows all of the virtual hoppings to 
add in parallel, thus overestimating the hybridization en¬ 
ergy gain, but because the mulitconfigurational character 
of the HS state is weak, the difference between DFT+U 
and DFT+DMFT is slight. 
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FIG. 4. Left panel: Energy splitting Ehs — Els as function 
of J for U = 5eV. Right upper panel; Energy splitting as 
function of U for J = 0.85eE. Right lower panel: Energy 
splitting as function of double counting parameter U' (Eq. 1) 
for U = 5eR and J = 0.85eR. Temperature is 387K in DMFT 
calculations. SDFT+U calculations were performed with pro¬ 
jector basis in VASP. 


Having established accurate energetics, we now turn 
to spectra and magnetic response. X-ray absorption 
spectroscopy experiments in which incident X-rays are 
tuned to the Fe L/// edge probe the empty d-states, re¬ 
vealing information about the electronic configuration of 
the Fe ions. Representative data^® are compared to our 
DFT-t-DMFT calculations in Fig. 2. The theoretical cal¬ 
culations (right panel) reveal that in the LS configuration 
the density of empty t 2 g states is very low; the eg spec¬ 
trum reveals a single main peak, with a small prepeak of 
t 2 g origin arising from a small probability of a LS d®. Al¬ 
ternatively, in the HS situation the two peaks correspond 
to transitions into the empty t 2 g and Cg states respec¬ 
tively. The calculated spectra (for example the t 2 g — eg 
crystal field splitting in the HS state) are in reasonably 
good agreement with the data although not all of the de¬ 
tailed structure away from the main peaks is reproduced 
(see Ref. 45 for a possible interpretation). 

Fig. 3 shows the susceptibility calculated from Eq. 2 
along with experimental results from Ref. 12. We see that 
the structure with long mean Fe-N bond length indeed 
has a Curie susceptibility y ~ 1/T permitting its identi¬ 
fication as the HS state, while the short bond (LS) state 
has a very small susceptibility. Experimental values are 
0.5 /xb and for LS and HS, respectively,^^ whereas we 
obtained around 0.004/iB for LS and 4.6~4.8 ^b for HS. 
Our calculations only include the d-electron contribution 
to x; the nearly quantitative agreement with experiment 
suggests this is the dominant contribution. 

Crucial to the DFT-I-DMFT formalism are the values 
of the interactions in the correlated subspace and the 
double counting correction. We have demonstrated a re¬ 
spectable degree of accuracy using the standard double¬ 
counting approach and accepted values for the on-site 
interactions. However, it is critical to understand the 


sensitivity of the results to these approximations. The 
left panel of Eig. 4 shows that the DET-I-DMET and 
DET-fU results for the LS-HS energy difference depend 
strongly on J, as expected since J is the term in the 
energy favoring locally high spin configurations. The 
magnitude of the slope is sufficiently large that a rela¬ 
tively small increase in the exchange from J = 0.85eE to 
J = 0.9eV would change the sign of the energy splitting, 
demonstrating the importance of precisely knowing the 
exchange. However it is significant that that generally 
accepted value J = 0.85eH yields an exchange splitting 
with the correct ~ O.leV order of magnitude. 

Eig. 4 also shows that three widely used SDFT-I-U 
methodologies yield a qualitatively different result as 
compared to DET-fU (and DFT-I-DMFT) in two key 
respects: the magnitude and sign of the slope and the 
value of the J = 0 intercept. The difference in J = 0 
values of the energy splitting shows that the SDFT func¬ 
tionals have an effective J built in to them; compar¬ 
ison to the J-dependent DFT-I-U results demonstrates 
that Jeff for SDFT(LSDA)-bU, SDFT(PW91)-bU, and 
SDFT(PBE)-|-U are approximately Jeff ~ 0.75eH, 
Jeff ~ 0.93eH, and Jeff ~ 1.05eH, respectively. The 
counterintuitive finding that the SDFT-I-U energy split¬ 
ting increases with increasing J can be traced to the spin- 
dependent double counting correction, which overcom¬ 
pensates the effect of the J in the interaction, increasing 
the energy splitting. This is reasonable behavior given 
that all of the spin-dependent exchange correlation func¬ 
tionals incorrectly predict the energy splitting to be neg¬ 
ative at J = 0. Therefore, careful analysis is required in 
the use of SDFT theories as a base on which to build a 
correlated calculation. 

The right upper panel of Fig. 4 presents the U depen¬ 
dence of the HS-LS energy splitting. We see that for each 
method, a U can be found that reproduces the measured 
energy difference, and the trends with U are similar in all 
methods, but the DFT-I-DMFT method gives the physi¬ 
cally correct splitting when a reasonable U is employed. 
Motivated by previous work on rare earth nickelates,^°’^^ 
we show in the right lower panel the effects of varying 
the double counting correction, setting U' ^ U in Eq. 1. 
For given U, J and U' , the DFT-I-DMFT procedure al¬ 
ways yields a smaller energy difference than the DFT-fU 
methodology. The dependences illustrated in Fig. 4 indi¬ 
cate that in order for the method to become truly predic¬ 
tive, improved theoretical understanding of the interac¬ 
tion parameters and double counting is required. Results 
presented here can serve as benchmarks for this endeavor. 

In summary, we have shown that DFT-I-DMFT with 
generally accepted interaction parameters produces en¬ 
ergetics, magnetic susceptibilities, and x-ray absorption 
spectra in reasonable agreement with experimental mea¬ 
surements on Fe(phen) 2 (NCS) 2 . The method involves 
a full self-consistency between the correlated subspace 
and the background, but the locality assumption basic 
to many solid-state applications of DMFT is here exact 
because there is only one correlated site. The ability 
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of DFT+DMFT to handle hybridization in a correlated 
environment is important for the success of the method. 
The ability to perform calculations for a range of temper¬ 
atures and structures revealed that electronic energy and 
entropy considerations do not account for the observed 
transition. The mean metal-to-ligand bond length is the 
key parameter controlling the spin state and the transi¬ 
tion is likely driven by phonon free energy considerations. 
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I. MAGNETIC SUSCEPTIBILITY 

Fig. SI compares the magnetic susceptibility of the 
large d{Fe — N) and short d{Fe — N) structures over 
a wide temperature range. The very small value and 
weak temperature dependence of x in the short d{Fe—N) 
structure demonstrates that this structure has a low-spin 
ground state, while the Curie-like temperature depen¬ 
dence of the susceptibility in the long d{Fe — N) struc¬ 
ture clearly shows that the local moment is well formed, 
even at the highest temperature studied here. 



d{Fe — N) structure an unphysically large J ~ 1.4eU is 
required to drive the transition. Thus, for a large window 
of J (0.5 eV to 1.5 eV) including physically reasonable 
values, DFT-I-DMFT is able to separate two spin states 
based on atomic structure. 


II. EE SITE ELECTRONIC CONFIGURATION 

Aspects of the Fe site electronic configuration can be 
directly measured in the DMFT calculations®^’®^. A 
measurements of the fraction of imaginary time that the 
impurity spends in a given configuration in the simula¬ 
tion yields the fractional weight of the configuration in 
the density matrix. Panels (a) and (b) of Fig. S4 show 
the occupancies of different Sz states. In the LS, 5*2 = 0 
is the most probable configuration, while for HS the most 
probable state is Sz = ±2.5. Orbital occupations (panels 
(c) and (d) in Fig. S4) also show the distinction between 
two states. In LS, t 2 g orbitals are close to full occupation 
and Bg orbitals are less than half occupied; but in HS, all 
orbitals are close to half occupation. 


FIG. SI. Temperature dependence of the Fe magnetic sus¬ 
ceptibility X from DFT-I-DMFT calculations in LS and HS 
structures with U = 5eF and J = 0.85eF. 

Fig. S2 shows the temperature dependence of the mag¬ 
netic susceptibility for a series of structures with different 
average Fe-to-N bond lengths. This figure demonstrates 
the existence of a critical value of the mean bond length 
Ri 2.10A above which the molecule is in the HS state and 
below which it is in the LS state. 

Fig. S3 shows the magnetic susceptibility calculated 
as a function of Hunds coupling J at a temperature T = 
387K slightly higher than room temperature for the short 
d{Fe — N) and long d{Fe — N) structures. We see that 
for very small J, the ground state of both the long-bond 
and the short-bond structure is LS. For both structures, 
as J is increased a LS-HS transition occurs. The critical 
value of J is r; 0.4eU for the long d{Fe — N) structure; 
this value is smaller than the physically reasonable J ^ 
0.9eU expected for Fe. On the other hand for the short 
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FIG. S2. Dependence of magnetic susceptibility on tempera¬ 
tures and average Fe-to-N bond length. 
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FIG. S3. Magnetic susceptibility as function of exchange J 
from DFT+DMFT calculations with U = 5.0ey and T = 
387K. 



FIG. S4. Spin (upper left and upper right) and (lower left 
and lower right) sector statistics for LS ((a) and (c)) and HS 
((b) and (d)) from DFT+DMFT calculations at t/ = S.Oel/, 
J = 0.85el/ and T = 387K. 


III. MEAN D OCCUPANCY AND DOUBLE 
COUNTING 

The LS and HS states differ in relative occupancies of 
the t 2 g and eg orbitals. This difference is the result of a 
competition between Hunds exchange J and crystal (lig¬ 
and) field splitting, which arises mainly from hybridiza¬ 
tion between Fe d and N-p orbitals. The crystal field 
splittings were calculated as the difference between av¬ 
erage energies of t 2 g and eg orbitals obtained from the 
on-site terms in the projection of the DFT Hamiltonian 
onto our MLWF basis. The calculated crystal field split¬ 
tings are shown in Fig. S5 and are clearly seen to be 
strongly influenced by the mean Fe-to-N distance 

Previous work®^^®® has shown that the occupancy Nd 



FIG. S5. Grystal field splitting as function of Fe-to-N dis¬ 
tance. 



FIG. S6. Dependence of effective magnetic moment pe// = 
obtained from on occupancy Nd of Fe d orbitals in high 
spin state at U = 5.0eU, J = 0.85eU and T = 387K. Differ¬ 
ent values of Nd are obtained by varying the coefficient U' in 
the double counting formula of the main text. 

of the correlated orbitals provides a useful characteri¬ 
zation of the beyond-DFT many-body state, and can 
be varied by changing the double counting coefficient 
U'®®. Nd is also related to the effective magnetic moment 
of the HS state and thus can be inferred from experi¬ 
ment. We have computed the effective magnetic moment 
Me// = by performing fully charge self consistent 

DFT+DMFT calculations using different values of U'. 
Measurements in the HS state indicate peff = 5/tB®^. 
As we can see from Fig. S6, for all Nd we studied in the 
HS state (4.7^5.5), our calculated effective magnetic mo¬ 
ments are smaller than the experimental value. The stan¬ 
dard FLL double counting U' = U yields an Nd ~ 5.25 
and an effective magnetic moment close to the maximum 
value and to experiment. 

We find that the d-occupancy is much higher in the 































FIG. S7. Nd from DFT+U and DFT+DMFT calculations performed for short d{Fe — N) (calculation restricted to LS state) 
and long d{Fe — N) (calculation restricted to HS state) structures, with various U, J and U’. Temperature is 387K in DMFT 
calculations. 



d(Fe-N) (A) 


FIG. S8. Total energy of LS and HS from DFT+U calcula¬ 
tions with and without charge self-consistency. U and J are 
fixed at 5.0 eV and 0.85 eV. 


it is to the interaction parameters. Of particular interest 
is the J dependence of Nd in the structure with the long 
d{Fe—N). In this structure, within the Hartree-Fock ap¬ 
proximation a level crossing between two HS states with 
different Nd occurs at a J ~ 0.5eV. A simplified picture 
of this level crossing is that as J is decreased it becomes 
favorable to add one electron to a minority-spin t 2 g level, 
thus increasing Nd and reducing the moment by Ri i. We 
emphasize that this metastable high Nd, high spin state 
is found in the HF calculations but not in the DMFT 
calculations. A truly low-spin state with Nd ^ 6 is also 
metastable within the Hartree-Fock approximation, but 
has higher energy than the state shown here. Level cross¬ 
ings in d-occupancy related to different spin states were 
also reported in a recent study of the heme molecule.®® 


LS state than in the HS state. The physics is that in 
the HS state the longer d{Fe — N) distance weakens the 
hybridization while the minority spin orbitals are pushed 
to a very high energy and thus hybridize much less with 
the ligands. Fig. S7 shows Nd computed with Hartree- 
Fock (DFT+U) for the short d{Fe — N) (“LS”) and long 
d{Fe — N) (“HS”) structures. Note that for the long 
d{Fe — N) structure the HS state is at least metastable 
within the Hartree-Fock approximation at all J which 
were studied, unlike the DMFT calculation where a HS- 
to-LS transition occurs as J is reduced below r; 0.4eU in 
the long d{Fe — N) structure. This is why the DMFT 
curves only extend down to J = 0.6 in the left panel. 

We see that where the HS state is stable in both meth¬ 
ods, DFT+U and DFT+DMFT give reasonably consis¬ 
tent pictures for the value of Nd in the two structures, 
and that Nd is much more sensitive to the structure than 


IV. THE EFFECT OF CHARGE 
SELF-CONSISTENCY 

Charge self-consistency can have a substantial impact 
on the total energy if there is a major rearrangement of 
charge from the initial density which was used to build 
the Kohn-Sham potential. The starting point of our 
DFT+DMFT calculation is a spin-unpolarized DFT cal¬ 
culation, which is inherently low spin. While this density 
will be very similar to the converged LS solution, it will 
be quite different from the converged HS solution. This 
is due both to the change in orbital occupancy and the 
overall large change of Nd across the LS-HS transition. 
In Fig. S8 we illustrate this effect by comparing the to¬ 
tal energy within DFT+U both with and without charge 
self-consistency. As expected, the changes are relatively 
small for the LS state and relatively large for the high 
spin state. 
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